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Abstract. Semiclassical approximations are implemented in the calculation of 
position and width of low energy resonances for radial barriers. The numerical 
integrations are delimited by r/r life <C 8, with r the period of a classical particle 
in the barrier trap and ruf e the resonance lifetime. These energies are used 
in the construction of 'haired' short range potentials as the supersymmetric 
partners of a given radial barrier. The new potentials could be useful in the 
study of the transient phenomena which give rise to the Moshinsky's diffraction 
in time. 

1 Introduction 

The study of resonant scattering processes has received considerable attention in contem- 
porary physics. Resonances are experimentally observed in atomic, nuclear, and particle 
physics so that diverse theoretical models have been proposed for their analysis over the 
years (see, e.g. [U(2])- In a simple picture, a resonance is a special case of scattering state 
for which the 'capture' of the incident wave produces delays in the scattered wave. The 
'time of capture' can be connected with the lifetime of a decaying system composed by the 
scatterer and the incident wave. Then, the resonance state is represented by a solution of 
the Schrodinger equation associated to a complex eigenvalue e and satisfying purely out- 
going conditions (Siegert functions) [3]. These functions are not finite at r — )• 00, so that 
they are not admissible as physical solutions into the mathematical structure of quantum 
mechanics. Some approaches extend the formalism of quantum theory to the wider back- 
ground where the Siegert functions are more than a convenient model to solve scattering 
equations [3HB]. Notwithstanding, the 'unphysical' behavior of the Siegert functions has 
been relevant in the construction of complex supersymmetric partners of a given poten- 
tial [7H9] (see also rjUJQT]). The main problem is to evaluate the complex point e up to a 
reasonable precision; the binding energy E = Re(e) and the lifetime r^e — ~l/(2I m (e)) 
of the decaying composite are then automatically determined. Below, the derivation of 
analytical expressions for E and THf e is discussed such that integrations are achievable up 
to the precision delimited by r/rufe <C 8, with r the period of a classical particle in a bar- 
rier trap. The results are used to get Darboux-deformations of radial barriers presenting 
'hair' over the top. A stronger resonant phenomenon is expected to be associated with 
these new interactions. 
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The avoiding of the asymptotic divergence of the Siegert functions is usually faced with 
r-complex coordinates r = pe tX , A > 0. The Siegert functions are then easily obtained by 
numerical integration, though they become eigenfunctions associated to complex eigen- 
values of a non-Hermitian Hamiltonian, as discussed in Sect. [2J To improve the numerical 
approximations it is feasible to include variations in the complex r-plane [12J, as is briefly 
outlined in Sect. 12.21 The latter approach, however, depends on one's ability to guess a 
trial function as a reasonable approximation to the actual wave function [13J. Variants of 
the WKB method are also available to get a more direct mechanism of integration [T3HT6]. 
Indeed, in most cases, bound and resonance energies correspond to poles of the S-matrix, 
since they are zeros of the Jost function involved. Such connection suggests that the 
bound energy techniques can be extended to the resonance case. Section [3] deals with 
the application of semiclassical approximations to calculate E and T in the one-channel, 
s-wave situation. The embedding of the energies into the complex e-plane due to small 
variations of the related wave-function is necessary to ensure the analytical continuation 
to the lower half-plane. As a test, the approach is applied to the simple scattering problem 
reported in [12J. Our results are in good agreement with those reported in e.g. [T2|[T7]. 
Applications to more elaborated potentials (as those reported in [16]) are straightforward. 
The implementation of simple and double complex Darboux transformations is discussed 
in Section HJ Interestingly, 'haired' barriers are found to be the supersymmetric partners 
of the conventional radial ones. Connections to the time delay problem in quantum me- 
chanics are transparent, so that the haired barriers could be of interest in the studying of 
the Moshinsky's diffraction in time. The paper is closed with some concluding remarks. 



2 Resonances and radial barriers 

Consider a spinless particle in a spherical, suitably short range potential U(r). In dimen- 
sionless form, the reduced radial Schrodinger equation is given by 



d 2 £(£+!) 



dr 2 r 2 



V(r) + k 2 



ip(r) = 



where R(r) = ip(r)/r is the radial part of the complete solution 0(f). When r — > oo, 
the effective potential vanishes and ip behaves like combinations of e ±lkr . Let us take 
k = Ke ia 6 C, a G R, k := \k\ and r > to write ikr = krcosa — Krsina. Henceforth 
the plane wave functions 

(p ± (r, k, a) = e ±ikr = ^""e^™ 

can be analyzed in terms of a. Three general cases are distinguishable: Scattering energies. 
If a = then k = n and the complex functions <p± = e ±lKr oscillate with finite amplitude 
at large distances. Bound energies. If a — 7r/2 one gets k = in, and the real function 
is divergent, while (p + = e~ Kr becomes zero for large values of r. Resonance energies. If 
a = —(3, with < < ^, then k is in the fourth quadrant of the complex /c-plane and the 
oscillation amplitude of the complex function <p_ (<£>+) decreases (increases) exponentially 
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as r — > +00. In any case, the solutions of (pQ) that actually behave like <p± are the Hankel 
functions hf(z), for which hf(z) — > e ±l ( z - i7r / 2 ) as 2; — > 00 (see, e.g., [13]). Therefore, 
solutions corresponding to bound or to resonance states are picked out from the linear 
combinations of hf(z) such that only outgoing waves exist. That is, these solutions satisfy 
the Siegert condition, 

lim _i-_^( r ) = ik, ke C. (2) 

Thereby both bounded and resonance functions behave as ip + at large distances. Clearly, 
these solutions must also be regular at the origin and smooth in between. 



2.1 Complex-scaling method revisited 

Let us take the transformation H — > \JHU~ 1 = Hg, with 9 a dimensionless parameter 
and the operator U such that r — > re l6 . Then ikr — > ikre %e and the bound wave-functions 
<f+(r, in, |) become exponential decreasing complex functions for large values of r provided 
that — ^ < 9 < |. In turn, the Siegert functions are transformed as follows 

y? + (r, k, -p) -> e lKrco s(e-/3) e -™(e-/3) ; < /3 < ^. (3) 

Thereby, <p+{r, k, —(3) is mapped into a bounded function if 9 — /3 > 0, i.e. for # e (0, |). 
As regards the scattering states, we have: 

To preserve the original form of scattering wave-functions the kinetic parameter k = k 
must be also modified, thus k — > ne~ ie . This last transformation induces a rotation of 
the positive real axis of energies in the clockwise direction by the angle 29 : E oc k 2 — > 
K 2 e~' l2e oc Ee~ l2e . That is, the rotated scattering energy is complex, e = €r — i^, with 
e# = E cos(26*) and ^ = ^8^(26*). In summary, the scaling operator [/ produces an 
embedding of the energies E G M. into the complex e-plane by a rotation of E > such 
that: 



(1) The bound state energies are preserved 

(2) The cut is rotated downward making an angle 29 with the real axis, and 

(3) The resonances are exposed by the cut. 

In general, the new 'complex eigenvalues' are independent, so that the resonance phe- 
nomenon is associated to the discrete part of the complex-scaled Hamiltonian Hg. Finally, 
considering the operators R and EE, with [R, II] = il, one realizes that U = e~ eRU is the 
scaling operator we are looking for. Indeed, since Hg = e~ 2l9 U 2 + V(e t9 R) is such that 
Hg 7^ Hg, the 'regularized' Siegert functions (J3J) are square-integrable eigenfunctions as- 
sociated to complex eigenvalues e of the non-Hermitian Hamiltonian Hg. 
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2.2 BBJS complex-coordinate analysis 



The repulsive potential 

V(r) = V r 2 e- Xr (4) 

represents a radial barrier of maximum height V max = 4Vo/(Ae) 2 at tq = 2A" 1 , which 
admits narrow shape resonances in low-energy scattering. Introduced in 1974 by R. A. 
Bain, J. N. Bardsley, B.R. Junker and C. V. Sukumar (BBJS), this potential is widely used 
in the testing of diverse approaches developed to calculate resonances (see e.g. [T2|IT7l-20j). 
In their paper, Bain and coworkers present a combination of the complex-scaling method 
with variational principles to analyze the resonance phenomenon. Indeed, they found a 
single resonance eo = 6.8722 — iO. 025549 (a factor 2 must be considered in [T2]) for the 
potential (J4]) with A = 1 and Vo = 15. To depict the BBJS procedure consider equation 
(PQ). After a complex scaling transformation, we get Mgip = e 2ie [Hg — k 2 }ip = 0. The next 
step is to assume that the integral 

POO 

1=1 jjj(r)M g ip(r)dr 
Jo 

is preserved for diverse values of 9 and small variations of ip(r). Introducing a trial 
function f e , parameterized by a set of variable numbers Ci,...,c n , the approximation 
f t (r;ci, . . . ,c n ) ~ ip(r) is better for small variations of / with respect to each of the 
parameters. That is, every one of the expressions j¥- ~ becomes an equality if f e = 
To state precisely: the BBJS method allows for the determination of approximate values 
of k, leading to points ?=fc 2 eC which are in the vicinity of the resonances e = E — i~ 
we are interested in. 



3 The WKB method for resonances 

In this section we formulate an extension of the WKB method to include the calculation 
of complex- valued energies e = E—i~. For simplicity, we shall focus on low energy s-wave 



shape resonances. We use the abbreviation p(r) = \/\T\, with T = E — V(r), the kinetic 
parameter for a given energy E. The WKB wave function in a classically accessible region 
(T > 0) and in the nonclassical domain (T < 0) is, respectively, written as follows 

V>(r) = -±=e**<*r\ ^(r) = -L=e ±n ^\ (5) 

\M r ) vw) 

where 

rr 

W(c,r) = / p(r) dr, f2(c, r) = iW(c,r). (6) 



If p(a) = p(b) = 0, the roots r = a and r = b are the classical turning points for energies 
below the barrier (see Figured]). The origin is a fixed turning point so that (0, a) and 
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(b, +00) are the classically allowed regions while (a, b) is the nonclassical domain. The 
connection formulae can be summarized as follows 



1 / . 7r \ smwe 

—j== cos W{c, r)-- + w) ^i) ->• == 

\/p{r) v 4 / v^v 



1 cos wt 

2 y / p(r) 



(7) 



where w is a real parameter defined by the turning point which is under analysis. We 
want to get the quantization rule for the resonance energy e = E - i\ by imposing the 
Siegert condition and looking for regular solutions at the origin. 




Figure 1: A typical short-range radial barrier {solid curve). For E given [dashed line), the 
classical turning points define two classically accessible regions (E > V), in the sequel labelled 
1 and 3 from left to right, and a nonclassical domain (E < V) labelled hereafter region 2. 



Let the solution in the nonclassical domain be written in the form 

-fi(r,6) 



(8) 



with A and B constants. The additivity of the Q-functions Q(a,b) = Q(a,r) + Q(r,b) 
leads to 

A n ,„ M tt B 



_ e -n(a,6) e w(r,6) _|_ ^_ 



-fi(r,6) 



This last solution is connected to that of region 3 via (JZj) with w — ~: 



2z5 

V?2 ->■ — P== I COS 



W(6,r) - 



7T' 



4J 45 



e Q(a ' b) sin 



7T 

4J 



Now, imposing AB = Ae n( - a ' b ^ we arrive at the plane wave 



A < r 

^exp\-Q(a,b) +1 W{b 



r 

; 4 



2 v / p(r 

In a similar manner, the connection to region 1 (w — 0) leads to 



.4 



cos 



7T 



W(r,a)-- 



V 2f ^ sin 
4 



7T 



W(r,a)-- 



(9) 



(10) 



(12) 



At the origin one must get <pi(0) = 0. However, this condition is not trivially fulfilled by 
( I12p so that we first constrain the W^-function to satisfy 



W(0,a) =[n + -)ir, 



(13) 



and then we impose the condition 



-2Q(a,b) 



< 1. 



(14) 



That is, we look for a function (fT2|) such that for small values of r the real part goes to zero 
faster than its imaginary counterpart. The imaginary part of ifi(r) is then expected to 
produce the embedding of the energy eigenvalues into the complex e-plane by the addition 
of a small (negative) imaginary term to each of the energies: 



r 



R3 E ^ £3 e = E + 5E, 5E = T > 0. 

This embedding is relevant since W is displaced to W + 5W, where 



SW(0,a) 



p 1 {r)dr 



T tT 

SE = -8E = -%—. 

4 8 



(15) 



Here r corresponds to the period of a classical particle moving harmonically from the 
origin to the turning point a. Thus, r is a quadrature depending on the constants of 
integration E and a. Assuming SW(0,a) <C 1, the variation of equation (j!2p shows that 
Sipi(r) is the correction for (pi(r) to be zero at the origin of the complex plane. In other 
words, (p(0) + 5tp(0) = provided that 



-2fi(a,6) 

5W{0,a)+i = 0. 

From ffl5l) we finally get an analytical expression to calculate the resonance width: 

Y e -2fi(a,6) 



(16) 



(17) 



The time spent by a trapped quantum in the tunnelling of the radial barrier is therefore 
given as follows: 

1 



dE 



Remark that our approach is useful for < rr C 8; otherwise ffl4l) is not satisfied. 
Notice also that this result is consistent with our assumption that SW(0,a) 1. Now, 
the turning points a and b move farther apart as r — > and, as a consequence, the factor 
e -2fi(a,&) deceases exponentially. Therefore T is narrower as a — > 0. On the other hand, 
a and b move closer as E — > V max , so that T increases as a — > b up to E ~ V m3X , where 
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BBJS 




KLM 




Extended WKB 




15[08.1201] 


6.8722 - i2.5549(-2) 


06.994 


-z2.787(- 


-2) 


07.01129 - z3.71716(- 


-2) 


30[16.2402] 




11.104 


-zl.321(- 


-4) 


11.05705 - il.41354(- 


-4) 


45 [24.3604] 




14.288 


-z6.840(- 


-7) 


14.21889 - i7.05499(- 


-7) 



Table 1: The single resonance eo for potential with A = 1 and the indicated values of 
Vo[V max ]- A factor 2 must be considered in the data originally reported by BBJS in [12] and 
KLM in [IT]. 



the simple WKB method does not apply [H]. For energies E > Vmax, the turning points 
diverge into complex values and the approach is not directly applicable (see, however, [16] 
and P2]). 

In summary, our method gives better results for energies which lay deep with respect to 
V m3iX ; improvements can be achieved by considering higher-order phase integral methods 
[2T] . Quite remarkably, the approach presented here could be applied in the study of time 
delay where transient effects are known to be relevant [22]. That is, the method could be 
useful in giving a quantum definition to the difference between the time to traverse the 
barrier and the time of going from a to b as a free particle. Particularly if the energy E 
is below Vmax, in which case, classically the particle could not arrive at the region 3 [23] 
(see also [24"J). 



3.1 Application: BBJS potential scattering 

We have derived two expressions to calculate resonances. Equation ( JT3l) corresponds to the 
WKB quantization and localizes the position E of energies fulfilling the Siegert condition. 
Equation (fT7|) . on the other hand, defines the width F of the resonance up to the precision 
established by < rT <C 8. We can go a step further in our approach by noticing that 
potential (J4J) admits a series of low-energy resonances, in correspondence with the strength 
Vo. For getting 'low-energy' resonances, we adopt the criterion of selecting those complex 
eigenvalues e whose real part is smaller than, or equal to Vmax- Thus, for a = r such that 
V(a) = Vmax, we want W in (1131) to have a maximum at a = r . This condition is that 
the integrand of PV(0,ro) involve the longest distance between E and V(r) for each point 
in [0,ro). The solution of this extremum problem enables us to identify the number n 
of resonances in terms of Vo, as requested. Indeed, the straightforward calculation shows 
that potential (@J admits n resonances, provided V is limited as follows: 

V (n-l)<Vo<Vo(n), nGN, (18) 

where 2 

%)=A|n + ! , A>=f^, n = 0,l,2,... (19) 
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and 

7= / \fl - z 2 e 2 ^- z )dz. (20) 
Jo 

A numerical integration gives 7 ~ 0.55098, so that Pq ~ 15.014 and potential fll]) admits 
a single resonance if 8.44539 < Vo < 45.9804. The strength Vq = 15 used in the BBJS 
paper [T2] is clearly in this category. The same potential is also analyzed by H.J. Korsch, 
H. Laurent and R. Mohlenkamp (KLM) in their study of the Milne's differential equation 
[IT] . Table [1] shows the good agreement of our results with those reported in |T2j and 
[T7] for a single resonance. Notice that eo is such that its real part is close to V m3iX = 
8.12012. Figure [2] shows the Siegert function obtained as a numerical solution of the 
corresponding Schrodinger equation for the extended WKB value of e . The closeness of 
the single BBJS resonance to V max is also observed for the 'highest' resonance in each of 
our studied cases. Namely, the highest resonance is such that Re(e n ) < Re(e]y)Wn < 
N and Re(eTv) ~ V max . For instance, Table |2] includes the values of the unique five 
resonances belonging to potential (TjJ for A = 1 and V = 350. The highest one is such that 
Re(e 4 ) ~ V max = 189.46939. For the lowest resonance, on the other hand, the bigger the 
difference V max — Re(e ), the narrower the width r , as was noticed in the previous section. 
According to Table El for example, Vo = 350 is such that V max — Re(eo) = 143.02639 and 
r /2 = 4.2321 x 10~ 36 . The strength V = 15 gives in turn V max - Re(e ) = 1.10883 and 
r /2 = 0.03717. Under these conditions the classical motion of a trapped particle takes 
place for larger periods r in the former than in the second case. 



-1 - 




5 10 15 20 



Figure 2: The real (solid curve) and imaginary (dotted curve) parts of the numerically integrated 
Siegert function belonging to the single resonance of V(r) = 15r 2 e~ r . 

As a final illustrative example of the method we include the anharmonic spherical 
oscillator V(r) = \kr 2 — gr N , k > 0, iV = 3, 5, . . . discussed by Mur and Popov in [16]. 
This potential has a barrier of height V max = ^p-k N/ (Ar_2) (gN) ~ 2/ (Ar_2) at r^ 2 = ±. 
Considering s- waves and the parameters N = 3, k = 800 = 2g (V max = 59.25925), the 
potential possesses a single low-energy resonance at e = 47.0105 — iO. 861937, which could 
be also tested by using the BBJS or the KLM methods. Hotter resonances for which 
E > Vmax can be analyzed with either the Mur-Popov approach or the KLM one. 
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Re(e) 


-Im(e) = T/2 




46.4430 


4.23210 x 10" 36 


ei 


95.9768 


2.91608 x 10~ 21 




135.8367 


1.23305 x 10~ n 




167.0870 


5.98751 x 10~ 5 


e 4 


188.5395 


0.71602 



Table 2: Extended WKB values of the unique five resonances for potential (j4|) with A = 1 and 
V = 350 (F max = 189.46939). 

4 Susy transformations and concluding remarks 

The Darboux transformation 

V{r) = V{r) + 2(5'{r) (21) 

is useful in many branches of the mathematical physics (see e.g. [25J ) . Of particular 
interest in quantum theories, this transformation supports the mathematical structure of 
the supersymmetric approach (for recent reviews, see [26H30] ). Moreover, the axioms for 
Hermitian operators and real spectra can be abolished in some situations (cf. [714TT |[3lT[35] 
and [36JI37]). Two potentials V and V are said to be supersymmetric (Susy) partners if 
the /9-function in (12ip is a nontrivial solution of the Riccati equation, 

- (3'{r) + f3 2 {r) = V(r) - e. (22) 

In such a case, the spectrum of V is the same as that of V, occasionally extended by 
an additional point e which could be complex. The non-linear equation ( )22l) is linearized 
to the Schrddinger equation Hip e = eip e by means of the logarithmic transformation 
(3(r) = — j- In ip e (r), with not necessarily in £ 2 (DomH). For tp e satisfying the Siegert 

condition (j2J), one easily verifies that potential V inherits the behavior of V at infinity. It 
is a custom to take -^-functions with at most a single zero, provided this root is one of 
the borders of Dom(H). In the present case, ip t is the Siegert function belonging to one 
of the previously derived resonance energies. Thereby, V is complex-valued and singular 



V 


n 




Re(e ) 


-Im(eo) = T /2 


15 


1 


8.12012 


7.01129 


0.03717 


60 


2 


32.48046 


16.8584 


4.86902 x 10~ 9 


150 


3 


81.20116 


28.8664 


2.14549 x 10~ 19 


250 


4 


135.33528 


38.5165 


1.97227 x 10~ 28 


350 


5 


189.46939 


46.4430 


4.23210 x 10~ 36 



Table 3: The extended WKB values of the 'lowest' resonance eo for different strengths Vq of 
potential ^ and A = 1. The number of resonances n is in correspondence with the condition 

GHD. 
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at the origin. A second Darboux transformation, using this time V and the complex 
conjugate of /3, gives rise to a new Susy partner of V, which is a real function. Also this 
new potential inherits the initial spectrum and is a regular function in DomH, as shown 
in Figure EJ 




Figure 3: Haired second order super symmetric partner of potential V(r) = 15r 2 e r . Each of 
the twice Darboux-distortions (couple of hairs) is a smooth curve, as shown at the right. 

The new radial potentials exhibit 'hair' along the negative slope, which induces stronger 
resonant phenomena. Since the Siegert function oscillates, while its amplitude increases 
exponentially (see Fig. [TJ, for large distances the amplitude of the /3-function decreases 
up to — ik. Hence, according to the Siegert condition (j2J), the complex double Darboux- 
distortions (couple of hairs) are cancelled as r — > oo. The same phenomenon is presented 
in square wells, where analytical expressions for E and V have been obtained and the 
number of hairs depends on the excitation of the resonance [HIES] ■ It is reasonable to as- 
sume that these haired potentials induce delays on the scattering states which are longer 
than the delay associated with their supersymmetric partners. It is then interesting to 
analyze the transient phenomena in the scattering process of these new potentials. In 
this way, the supersymmetric quantum mechanics could be connected to the Moshinsky's 
diffraction in time through these haired potentials. 
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